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Abstract 

rr The heat transfer coefficient between a molten charge and its surroundings in a Bridgman 
furnace was determined using a novel approach utilizing in-sitv temperature measurement. The 
ampoule containing an isothermal melt was suddenly moved from a higher temperature zone 
to a lower temperature zone. The temperature-time history was used in a lumped-capacity 
cooling model to evaluate the heat transfer coefficient between the charge and the furnace. The 
experimentally determined heat transfer coefficient was of the same order of magnitude as the 
value estimated by standard heat transfer calculations. 

1 Introduction 

A variety of directional solidification techniques are used for preparation of materials, espe- 
cially for growth of single crystals. These include the vertical Bridgman-Stockbarger, horizontal 
Bridgman, gradient freeze, and zone melting techniques. It is time consuming and costly to experi- 
mentally determine the optimal thermal conditions of a furnace utilized to grow a specific material. 
Hence, it is desirable to employ analytical and numerical models to assist in determining the optimal 
thermal conditions for a specific growth system, e.g. Bridgman technique [1-6]. Such calculations 
are handicapped by limited knowledge of the growth environment’s thermal characteristics. 

The heat transfer coefficient, defined as the ratio of the heat flux to the temperature difference 
between the material and the furnace, is an important thermal parameter in a growth system. The 
heat transfer coefficient manifests itself in the heat transfer models as the Biot number hR/k, where 
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h is the heat transfer coefficient, R is the sample radius, and k is the thermal conductivity of the 
sample. The Biot number may be regarded as the ratio of the ease of heat exchange with the 
furnace to heat conduction through the charge. 

Chang and Wilcox [1] showed that increasing the Biot number in Bridgman growth affects the 
position and shape of the isotherms in the furnace. The sensitivity of interface position to heater 
and cooler temperatures is greater for small Biot number, i.e. small ampoule diameter, ineffective 
heat transfer from the ampoule surface (small value of h), and high thermal conductivity of the 
material. 

Fu and Wilcox [2] showed that decreasing the Biot numbers in the heater and cooler of a vertical 
Bridgman- Stockbarger system results in isotherms becoming less curved. The planar isotherms lie 
in the lower portion of the adiabatic zone when the heater’s Biot number is larger than the cooler s 
Biot number. Increasing the cooler’s Biot number moved the position of the planar interface toward 
the upper section of the adiabatic zone. 

Although the heat transfer coefficient may be estimated from heat transfer principles [5,7], 
there are considerable uncertainties making an experimental value preferred. Here, we report 
a novel approach to experimentally determine the average heat transfer coefficient between the 
growth material and the furnace. This is accomplished by in-$itu temperature measurement of 
a transiently cooled object, i.e. melt or solid contained in an ampoule. In this technique, an 
isothermal charge held at temperature T 0 is suddenly moved to a chamber at temperature T^. 
The temperature in the charge as a function of time is measured using a thermocouple and the 
data obtained are used to calculate the average heat transfer coefficient between the molten charge 
and its surroundings using a lumped-capacity model. A GaSb charge was used to demonstrate the 
method. The validity of the lumped-capacity model was examined by solving the one-dimensional 
transient heat transfer problem between a rod and its surroundings and comparing the results 
of this analysis with the lumped-capacity solution. The one-dimensional transient problem takes 
radial temperature gradients in the rod into account. 

2 Lumped-Capacity Model 

The value of the heat transfer coefficient depends on the geometry of the system as well as the 
physical properties and temperature of the material and its environment. A simple but important 


method based on a lumped-capacity solution [8] can be utilized to determine the average heat 
transfer coefficient between an object and its surroundings from the transient cooling of the object. 
This analysis assumes that the object is isothermal. This object at temperature T 0 is introduced 
suddenly into an environment at temperature To©. The heat transfer coefficient is calculated from 
the change in temperature of the object as a function of time. 

An analytical temperature-time relationship of the cooling process is obtained by equating the 
change in the total enthalpy of the object to the net heat flow from the object to the surroundings 
during time interval dt: 

d [cpV(T - Too)] = -tA(T - Too)dt, (1) 

where the symbols are defined in the table of nomenclature. Assuming c, p, and h to be constant, 
integration of equation (1) from T 0 to T over time t gives: 

hAt 

In (T - Too) = In (T c - T^) - (2) 

Transforming the above equation into a non-dimensional form yields: 

In# = (— 2Bi)r, (3) 


where Bi=hR/k, 6 =(T-T 00 )/(T 0 -T 00 ), and r = at/R 2 . Thus, the Biot number is found from the 
slope of a In 6 versus r plot and the heat transfer coefficient is given by: 



( 4 ) 


In our experiments, the object consisted of a molten GaSb charge contained in a quartz am- 
poule (discussed in detail in the experimental section). The environment was a vertical Bridgman- 
Stockbarger furnace. The ampoule wall added a resistance to the heat transfer between the charge 
and the furnace. To include the effect of such resistance on the heat transfer coefficient, some mod- 
ifications of the thermophysical properties in equations (1-4) were necessary. (A similar approach 
was undertaken by Naumann [5].) A mass weighted effective thermal conductivity and thermal 
diffusivity are defined as follows: 
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where all parameters are defined in the table of nomenclature. The effective heat transfer coefficient 
then becomes: 


h«ff = 


Bi e ffk e ff 
R. ’ 


(7) 


where Ra is the outer radius of the ampoule and Biefr is found from the slope of the In 6 versus 
r e ff = <* e fft/Rl plot. 


3 Validity of the Lumped-Capacity Model 


In order to determine the validity of the lumped-capacity model for different ranges of the 
Biot number, a similar model was solved which takes radial temperature gradients in the sample 
into account. This model assumes that the heat flow is axisymmetric, the temperature in the rod 
is uniform in the axial direction, and the density, specific heat, and thermal conductivity of the 
rod are independent of temperature. It is valid for any values of the Biot number. The solution to 
this problem for the dimensionless temperature at the centerline of the rod is given in Carslaw and 


Jaeger [10] as: 


where the eigenvalues 


© _ T ~ Tco _ 2 f — exp 

° “ To - Too ~ 2 ^ Jo(r n )(r2 + Bi 2 ) p 

r n are roots of the equation: 



( 8 ) 


r n Ji(r n ) = Bi j 0 (r n ), (9) 

and Jo and J x are Bessel functions of order zero and one, respectively. For the range of Biot 
numbers studied, ten terms of the infinite series were determined to yield an accuracy of more than 
5 significant digits in 6. 


When the series in equation (8) is truncated after the first term, the following equation is 
obtained: 

In 0 X = In I (pz + Bi z j } " r ? T - (10) 

In order to determine when the series could be truncated after the first term, the error in trunca- 
tion had to be determined as a function of dimensionless time r and Biot number. The relative 
truncation error was defined to be (&i - 0)/6, where 6 is the exact solution which was taken as the 
value of 6 when the series was truncated after ten terms. The dimensionless time it past which 
the truncation error is less than 0.01 was plotted in Figure 1 for values of the Biot number between 
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0 and 0.9. For values of the Biot number between about 0.05 and 0.9, the critical dimensionless 
time is given by the b th order polynomial fit: 

Tent = 3.02 Bi 5 - 8.54 Bi< + 9.32 Bi 3 - 5.03 Bi 2 + 1.48 Bi - 0.0302. (11) 


This means that equation (10) is valid for all Biot numbers as long as the dimensionless time is 
greater than A plot of In# versus r should be linear after time Tan. 


If a linear regression analysis is carried out on experimental data at dimensionless times greater 


than T cr it, Ti is given by the square root of the negative of the slope. This value can then be used 
to calculate the Biot number by two different methods. A value for the Biot number can be found 
by solving equation (9). Also, the Biot number can be calculated from the intercept b 0 by finding 
the root of the equation: 


f(Bi) = e b ’ 


2Bi 


J.(ro (r 2 + bp) 


= o. 


( 12 ) 


The value of the Biot number calculated from the intercept is much more susceptible to 
experimental error. When doing a linear regression analysis, the value of the slope is much more 
statistically significant than the value of the intercept. This is because the intercept is at the edge 
of the experimental data. Also, the Biot number is much more sensitive to small errors in the 
intercept than in the slope. Figure 2 shows the dependence of the Biot number on the intercept b 0 
and slope m of equation (10). The solid line is a linear fit of the Biot number versus the intercept, 
and is given by: 

Bi = 4.24 b 0 . (13) 

The dashed line is a linear fit of the Biot number versus the negative of the slope, and is given by: 


Bi = 0.503 (-m). 


(14) 


Thus, the Biot number is about eight times more sensitive to experimental errors in the intercept 
than the slope. Another uncertainty when calculating the Biot number from the intercept is that 
the experimentally determined value of the slope must be used in equation (12) to find the intercept. 
This is an additional source of error. 


Equation (10) is similar to the lumped-capacity solution, except that the intercept is not zero 
and the slope is ~r? as compared to -2Bi from the lumped-capacity model. The validity of the 
lumped-capacity solution can be assessed by calculating the relative error caused by assuming the 
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Biot number equals the negative of the slope of a ln0 versus r plot divided by two. The relative 
error can be defined as (Bi-r?/2)/Bi, where Bi is the actual Biot number and Tj is calculated from 
equation (9). The relative error is plotted in Figure 3 as a function of Biot number and is given by 
the 2 nd order polynomial fit: 

Relative Error = -0.0373 Bi 2 + 0.248 Bi, (15) 

for Biot numbers between 0 and 0.9. The relative error is almost a linear function of Biot number 
in this region. 

This analysis suggests that care must be taken when applying the lumped-capacity model 
to experimental data. A linear regression analysis should only be done for values of r > t cnt- 
When the Biot number is large, most of the temperature change could occur before t and the 
error in assuming the slope equals -2-Bi becomes large. Furthermore, the experimental errors in 
temperature measurement become large when T approaches Too- This is best demonstrated by 
example. Let us consider the case when T e = 800°C, Too = 700°C, and the error in measurement 
of the temperature is 1°C. If the actual temperature is 750°C, the error in the resulting value of 
In# is only 3%. However, when the actual temperature is 701°C, the error is 15%. 

Another important point when applying the lumped-capacity model is that the straight line 
should not be forced through the origin. A two-parameter regression analysis should be performed, 
as suggested by equation (10). 

4 Experimental 

The lumped-capacity technique was demonstrated using an ampoule and thermocouple ar- 
rangement containing molten GaSb situated in the heater of a Bridgman-Stockbarger furnace at 
temperature T 0 . The ampoule was moved suddenly to a region of different temperature Too. The 
temperature versus time data were collected as the melt equilibrated to the new temperature and 
the lumped-capacity model was used to determine the heat transfer coefficient between the ampoule 
and furnace. 

A schematic diagram of the experimental apparatus is shown in Figure 4. It consisted of 
a 3-zone vertical Bridgman-Stockbarger furnace. The heating zones of the furnace were made of 
Kanthal heating elements embedded in Fibrothal insulation. The 5 cm long adiabatic zone was 
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fabricated from zirconia insulation. The zones’ temperatures were controlled to wi thing ±1°C using 
digital PID temperature controllers and SCRs. Quartz tubing was used as a liner in the furnace. 
Two type K thermocouples, inserted halfway into the heaters between the furnace wall and the 
liner, were used for control. Both ends of the furnace were plugged to eliminate the chimney effect. 

The 0.9 cm ID and 1.1 cm OD quartz growth ampoule, shown in Figure 5, was loaded with 
a 7 cm long GaSb charge, compounded from six-9s purity Ga and Sb in a rocking furnace for 5 
hours at 820°C. The temperature in the melt was measured using a 0.041 cm diameter grounded 
K-type thermocouple with a 310 stainless steel sheath and MgO as insulation (made by General 
Measurements). The tip of the thermocouple was positioned 3 cm into the ampoule at the center 
of the charge. 

The molten charge was allowed to reach thermal equilibrium with the surroundings prior to 
each experiment. For experiments HI and H2, the loaded ampoule was positioned in the lower zone, 
set at 890°C. The upper zone’s setting was at 800°C. After the thermocouple reading became stable, 
the ampoule was suddenly moved to the upper zone and held firmly. Meanwhile, the thermocouple 
output was collected using a data acquisition system. In experiment H3, the upper and lower 
zones were set at 870°C and 780°C, respectively. The ampoule in experiment H3 was initially 
positioned in the upper zone. After the temperature became stable, the ampoule was suddenly 
moved downward into the lower zone and held in the central region of this zone. The temperature 
was recorded versus time using a data acquisition system. 

5 Results 

5.1 Experimental Determination of Heat Transfer Coefficients 

Figures 6 and 7 show the actual thermocouple readings collected from the molten GaSb during 
experiments HI, H2, and H3. The logarithm of the dimensionless temperature 6 was plotted versus 
dimensionless time r e jf in Figures 8-10. A linear regression analysis using NCSS Statistical Software 
was performed to determine the slopes and intercepts of these plots. 

Figure 7 depicts the experimental data and the resulting linear fit for experiment Hi, where 
the GaSb charge was raised from the lower zone at 890°C to the upper zone at 800°C. The value 
of T 0 was 871.3°C and T^ was 792.3°C. The linear regression analysis was performed on the data 
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for r e ff between 0.1 and 18. There are two reasons that these limits for r e ff were chosen. The value 
of 7crit» beyond which the relationship between In# and r e g is linear, is less than 0.1. Also, the 
experimental data tend to bend upward beyond T e fl = 18. The reason for this is the magnification 
of experimental error at small values of 6. (Both of these explanations are discussed fully in 
section 3.) The resulting linear fit is 6 = (-0.142 ± 0.0006)r + (0.0697 ± 0.0059). The uncertainty 
values are 95% confidence limits on the slope and intercept, which can be transformed into 95% 
confidence limits on the Biot number. The lumped-capacity solution leads to a Biot number of 

0. 0712 ± 0.0003 determined from the slope. The two values of the Biot number determined using 
the model accounting for radial temperature gradients are 0.0725 ± 0.0003 from equations (9 k 10) 
and 0.0767 ± 0.0008 from equations (10 k 12). 

The same analysis as above was performed for experiment H2, which was a repeat of experi- 
ment HI described above. However, T 0 was 876.3°C and Too was 793.6°C in this experiment. The 
reason for the difference in these values from those of experiment HI was probably a difference 
in initial and final positions of the charge for the two experiments. The linear regression analysis 
was performed on the data for T e fj between 0.1 and 20 and the result is shown in Figure 8. The 
linear fit is 6 = (-0.138 ± 0.0004)r + (0.0407 ± 0.0042). The lumped-capacity solution leads to a 
Biot number of 0.0689 ± 0.0002 determined from the slope. The two values of the Biot number 
determined using the model accounting for radial temperature gradients are 0.0701 i 0.0002 from 
equations (9 k 10) and 0.0719 i 0.0005 from equations (10 k 12). 

The analysis was performed for experiment H3, where the GaSb charge was lowered from 
the upper zone at 870°C to the lower zone at 780°C. The value of T e was 854.5°C and Too was 
762. 5°C. The linear regression analysis was performed on the data for r e fj between 0.1 and 15 and 
the result is shown in Figure 10. The linear fit is 6 — (—0.171 db 0 . 0008 )t + (0.0522 ± 0.0073). The 
lumped-capacity solution leads to a Biot number of 0.0855 ± 0.0004 determined from the slope. 
The two values of the Biot number determined using the model accounting for radial temperature 
gradients are 0.0873 ± 0.0004 from equations (9 k 10) and 0.0903 ± 0.0012 from equations (10 k 
12 ). 

The heat transfer coefficients obtained for experiments HI, H2, and H3 are presented in Table 

1. There are three different values for each experiment. One is from the lumped-capacity model 
and two are from the model which accounts for radial temperature gradients. 
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5.2 Theoretical Estimation of Heat Transfer Coefficients 

The heat transfer coefficient between the ampoule wall and the furnace could also be estimated 
by a simple heat transfer formulation. The heat transfer coefficient is derived by summing the heat 
transfer by radiation and conduction through the air gap between the ampoule and the furnace 
wall in a concentric cylindrical system and equating it to the heat flux Q through the ampoule 
containing the growth materials: 

Q = E(T h -T.)-. ( F(Ti-T 1)+ £S^ (.6) 

The heat transfer coefficient determined from the above equation is: 

h = <reF(T 2 + T 2 )(T h + T a ) + r- 7 ^-. (17) 

R a ln^ 

The contribution of natural convection between the ampoule and the furnace was determined 
by computing the Grashof number in the air gap between the ampoule and the quartz liner: 

Q r _ g^»jr(Th — T»)(R] — Ra) Pair 2 ^g) 

n 2 

We estimated the Grashof number for our experiments to be 246. Since this is less than 1000 [11], 
we concluded that the contribution of natural convection to the heat transfer was not significant. 
Therefore, the convective heat transfer term was not included in equation (16). 

The view factor F = 0.9 was approximated using the view factor for finite-length concentric 
cylinders (i.e. furnace-ampoule) [13]. The ampoule temperature T a was assumed to be the same as 
the temperature of the molten charge, since the ampoule wall was relatively small compared to the 
charge diameter. The steady-state temperature reading in the charge Too and the temperature of 
the control thermocouple were used for T a and T^, respectively. The values of T a = 792. 5°C and 
T h = 800°C for experiments HI and H2, and Ta=762°C and Th=780°C for experiment H3 were 
used in calculation of the heat transfer coefficients. 

Using equation (17), the values of h calculated for experiments HI, H2, and H3 were 0.0220 
(W/cm 2 -K), 0.0220 (W/cm 2 -K), and 0.0189 (W/cm 2 -K), respectively. The estimated heat transfer 
coefficients are compared with the experimentally determined values in Table 1 . 
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6 Discussion 


Figure 11 depicts the comparison between the Biot numbers calculated from the lumped- 
capacity model (Bi = —m/2) and from the model including radial temperature gradients (Fj = v— ; m 
Bi = Tj Jj(Fi)/J 0 (ri)). The horizontal and vertical lines in the figure represent 95% confidence lim- 
its on the values of the Biot number. The three data points would fall along the solid diagonal 
line if the results from the lumped-capacity model agreed exactly with the theory which considers 
radial temperature gradients in the charge. The dotted line represents the calculated relative error 
between the two models and is given by equation (15). 

The Biot numbers calculated from the model accounting for radial temperature gradients were 
greater than those calculated from the lumped-capacity model. The fact that the three data points 
fall on the dotted line proves the relationship between the lumped-capacity Biot number and the 
radial temperature gradient Biot number is satisfied exactly. If the Biot number is determined 
using the simple lumped-capacity model, the correct Biot number can be calculated by finding the 
root of the rearranged form of equation (15): 

f(Bi) = 0.0373 Bi 3 - 0.248 Bi 2 + Bi - Bi L = 0, (19) 

for values of the Biot number between 0 and 0.9, where BiL is the Biot number calculated using 
the lumped-capacity model. 

In theory, the value of the Biot number obtained from the slope and equation (9) should be 
the same as the value obtained using the intercept. The solution to the problem which accounted 
for radial temperature gradients was for a charge made of a single material, while the charge in 
our experiments consisted of a molten semiconductor contained in a quartz ampoule. When radial 
temperature gradients were considered, the procedure of using mass weighted thermal properties 
probably led to underestimation of the radial temperature gradients present. The ampoule, which 
caused a discontinuity in the radial temperature profile, acted as an additional barrier to heat trans- 
fer that could not be properly accounted for by adjusting the physical properties. The difference 
between the values of the Biot number obtained from the slope using equation (9) and from the 
intercept represents a failure of the exact model. 

Experiments HI and H2 were almost identical, yet the Biot numbers calculated from each 
differ by about 4%. This difference could have been caused by a change in the position of the 
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charge within the furnace. The length of the heated zones were 15 cm, not overly long compared to 
the charge length of 7 cm. It is possible that one or both experiments were influenced by the charge 
being positioned in a region where the temperature changed with height in the furnace. This would 
cause axial temperature gradients in the melt which would violate the assumptions of the models 
used. However, differences in Biot numbers less than 10% are not significant when considering the 
fact that heat transfer coefficients depend continuously on position in a crystal growth furnace, and 
this technique is being used to get an overall value for the heat transfer coefficient in a zone. 

The Biot number calculated from experiment H3 was about 20% higher than those calculated 
from experiments HI and H2. Experiments HI and H2 measured the heat transfer coefficient in the 
upper zone, while experiment H3 measured that in the lower zone. The difference in heat transfer 
characteristics between the upper and lower zones was unexpected. These zones are constructed 
identically. The temperature rolls off at the bottom of the lower zone. It is possible that this 
temperature gradient led to axial heat transfer which increased the heat transfer coefficient. 

It is also possible that natural convection was present in the liquid during the experiments. 
However, the effect of natural convection would not be expected to be large in a low Prandtl number 
fluid like a semiconductor melt. 

The values of the heat transfer coefficient calculated experimentally differ by only 15% from 
simple theoretical estimations. This difference is quite small and could be due to inaccuracies in 
the values of emissivity and view factor. 

7 Conclusions 

A novel and practical approach was developed to determine the average heat transfer coefficient 
between a charge and directional solidification furnace. An ampoule containing molten GaSb and a 
thermocouple was suddenly moved from a high temperature region to a lower temperature region. 
The temperature decay versus time was recorded and used in a lumped-capacity model to determine 
the effective Biot number between the charge and furnace. The validity of the lumped-capacity 
method was assessed using a model which accounted for radial temperature gradients in the charge. 

It was determined that the lumped -capacity model is accurate within a relative error given 
by equation (15) for Biot numbers between 0 and 0.9. However, this is only if the linear regression 
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analysis on the In# versus r plot is carried out for values of time greater than a critical dimensionless 
time. This dimensionless time Tent is given by equation (11) for Biot numbers between 0.05 and 0.9. 
The straight line fit should not be forced through the origin, since the theory which accounts for 
radial temperature gradients predicts a slope and a non-zero intercept. For values of r greater than 
Tent, the Biot number can be determined exactly by obtaining T\ from the slope of a ln0 versus 
t plot. This value can then be used to calculate the Biot number from equation (9) or used in 
conjunction with the intercept of the In# versus r plot to find the Biot number from equation (12). 
However, the value of the Biot number determined using the intercept is much more susceptible to 
experimental errors. 

The difference between the experimentally determined Biot numbers calculated from the 
lumped-capacity model and the model accounting for radial temperature gradients was exactly 
as predicted by theory. If the Biot number is determined using the simple lumped-capacity model, 
the correct Biot number can be calculated by finding the root of equation (19) for values of the 
Biot number between 0 and 0.9. However, with the computing power and canned software available 
today, it might be easier to find the Biot number directly from equations (9 & 10). The Fortran 
functions DBSJ0 and DBSJ1 from the IMSL Library [19] were used by the authors to calculate the 
Bessel functions J 0 and Ji, respectively. 

The heat transfer coefficient between the ampoule wall and the furnace was estimated by a 

simple heat transfer formulation accounting for conduction and radiation across the air gap between 

the furnace and ampoule. Remarkably, the results differ by only 15% from the values calculated 

by the lumped-capacity method. 

% 
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Nomenclature 


A 

b„ 

Bi 

Bieff 

BiL 

c C 


C e ff 


F 

g 

h 


b e ff 


Jn 

k c 

k«ff 


kqa 

L 

m 

R. 

Rc 

Rf 

Ri 

t 

T 


Surface area of the charge (17.28 cm 2 ). 

Intercept of equation (10). 

Biot number *jjk. 

Effective Biot number between the charge+ ampoule and furnace 
Biot number calculated by the lumped-capacity method. 

Specific heat capacity of the molten GaSb charge calculated using 
thermal conductivity, thermal diffusivity, and density data (for GaSb 0.328 J/g-K). 
Mass weighted average effective specific heat capacity of charge+ampoule combination. 
c e ff = - y y ( 0.479 J/g-K for present experiments). 

Specific heat capacity of ampoule (c qa =l-19 J/g-K at 1075 K for fused silica [12]). 

View factor from furnace wall to the ampoule (0.9 for present configuration). 
Gravitational acceleration (980 cm/s 2 ). 

Average heat transfer coefficient between the growth material and 
the furnace (W/cm 2 -K). 

Effective average heat transfer coefficient between the charge+ampoule and 
the furnace (W/cm 2 -K). 

Bessel function of order n. 

Thermal conductivity of the charge (0.171 W/cm*K for molten GaSb [13]). 

Mass weighted effective thermal conductivity of charge+ampoule 
k.// = < 0146 W/cm K). 

Thermal conductivity of air (5.2xl0 -4 W/cm-K at 1100 K [14]). 

Thermal conductivity of ampoule wall (2.42xl0 — 5 T + 4.48xl0 -3 W/cm-K for fused 
silica [15], = 0.0308 W/cm-K at T=1075 K). 

Length of the charge (7 cm). 

Slope of linear equation. 

Outside radius of ampoule (0.55 cm). 

The charge radius (0.45 cm). 

Inside radius of furnace (3.81 cm). 

Radius of the furnace liner (1.0 cm). 

Time (s). 

Temperature at time t (°C). 
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I 


T, 

T. 

V v 

V c 

<*c 

<* e ff 

Pair 

e 

6 

V 

Pair 

Pav 

Pc 

Pqa 

(7 

r 

r n 


Steady state temperature of the melt (°C). 

Furnace temperature (°C). 

Ampoule temperature (°C). 

Volume of the section of ampoule containing the charge; 

Vqa = 7rL(R a 2 - Re 2 ) (2.2 cm 3 ). 

Volume of the charge (4.46 cm 3 ). 

Thermal diffusivity of charge (0.087 cm 2 /s for molten GaSb [13]). 
Effective thermal diffusivity of charge+ampoule (0.0628 cm 2 /s). 

Thermal expansion coefficient of air (0.00367 1/K [15]). 

Emissivity of the furnace (for Kanthal [18] € = 0.75). 

Dimensionless temperature ((T T^ )/(To — Too)). 

Viscosity of air (1.84xl0 -4 g/cm-s [15]). 

Density of air (1.19xl0~ 3 (g/cm 3 [15]). 

Average density of charge+ampoule, p av = (4-85 g/cm 3 ). 

Density of charge (5.98 g/cm 3 for molten GaSb at 800°C [17]). 

Density of ampoule wall, (2.28xlO -4 T + 2.273 g/cm 3 for fused silica [16] 
p qa = 2.586 g/cm 3 at 1075 K). 

Stephan-Boltzman constant (5.67xl0 -12 W/cm 2 -K 4 ). 

Dimensionless time (gr). 

Effective dimensionless time (^r)- 
n th eigenvalue. 
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Figure 2. The dependence of the Biot number on the intercept b 0 and slope m of equation (10). 
The solid line is the linear fit of the Biot number versus the intercept given by equation (13) and the 
dashed line is the linear fit of the Biot number versus the negative of the slope given by equation 
(14). 




Relative Error 



Figure 3. Relative error caused by assuming the Biot number equals -0.5 times the slope of a In# 
versus t plot. The line represents the 2^ order polynomial fit given by equation (15). 




Figure 4. Schematic diagram of vertical Bridgman-Stockbarger furnace. 





Figure 5. Schematic diagram of ampoule and 
system for in- situ temperature measurements. 
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Figure 8. Experimental data and the resulting linear fit (R 2 = 0.9993) for experiment HI. 
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Figure 10. Experimental data and the resulting linear fit (R? — 0.9991) for experiment H3* 




Experiment 

h (W/cm 2 -K) 

Lumped-Capacity 

Model 

Radial Model 

Theoretical 

Estimation 

Equations (9 & 10) 

Equations (10 & 12) 


0.0189 

0.0193 

0.0204 

0.0220 


0.0183 

0.0186 

0.0191 

0.0220 

H3 

0.0227 

0.0232 

0.0240 

0.0189 


Table 1. Comparison of the experimentally determined and theoretically estimated values of the 
heat transfer coefficient for experiments HI, H2, and H3. 




Biot from Equation (3) 



Figure 11. Comparison between the Biot numbers calculated from the slope using equation (3) 
(lumped-capacity model) and using equation (9). The horizontal and vertical lines in the figure 
represent 95% confidence limits on the values of the Biot number. The dotted line represents the 
calculated relative error between the two models and is given by equation (15). 



